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Novel solutions of Minkowski QED2+1 and QCD Schwinger-Dyson equations will be presented 
and discussed. The resultant propagators of confined degrees of freedom will be shown. 
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I. MOTIVATION AND INTRODUCTION 

bv 

In perturbation theory these are always the free propagators which enter into the Feynman diagrams hence the 
well known calculation trick known as the Wick rotation is always feasible. In this way the Euclidean and real world 
-Lorentz (Minkowski)- space results are related through a well defined analytical continuation. As a consequence of 
the perturbative technique, looking in momentum space, the Green's functions (GFs) of given quantum field theory 
have a real branch points whose positions can be traced by analyzing the matching propagator poles in their possible 
convolutions inside the Feynman integrals. The location of such branch points corresponds with physics - with opening 
' particle threshold production. For the momenta above threshold the GFs become complex with real and imaginary 
parts uniquely related by well known dispersion and unitarity relations. Assuming particles can freely propagate, the 
<~~~i appropriate propagators must have a real pole and the unitarity relations remain to be valid for S-matrix elements. 

This, of course completely follows the idea of LSZ reduction formula which relates S-matrix and Green's functions of 
I ■ given theory. 

Q-T Apart of genuine perturbation theory success in quantum electrodynamics and electroweak sector of the Standard 
Model, hadrons built of "light" quarks gain their masses through the approximate dynamical chiral symmetry breaking 
(DCSB) while there is simultaneously a nonperturbative mechanism which forbids colorless hadrons to disintegrate 
to any free quarks (or colored states in general). Similarly, the brehmstrahlung and consequential propagation of 
an on-shcll gluon could be similarly impossible as well. Clearly the phenomena of DCSB and color confinement lie 
beyond the access of the perturbation theory, which is based on particle-field duality and which starts to work with 
free propagators from the very beginning. 

On the other hand the GFs of confined fields are basic element to build hadronic wave functions form the first 
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. principle, i.e. using the original QCD Lagrangian degrees of freedom- the quark and gluon fields. The spectra of all 
hadrons have timelike signature P 2 > and the masses and and decays are identified from measurements performed 
in the real world- one time and three space coordinates forming the Minkowski space. The intentions of my talk is to 
' mention various non-perturbative methods which can be or have been used to evaluate GFs in strong coupling theories 
in Minkowski space. The natural framework is based on the setup of QCD Schwinger-Dyson equations [H, [1, El • 

The standard philosophy of nonperturbatively solved SDEs is common with the lattice theory. This is to start 
with Euclidean space definitions of Greens functions and analytically continue the Euclidean results to the Minkowski 
space. This is clearly always possible, even if the observed singularities do not reflect usual assumptions imposed by 
Wick rotation (4j. The SDEs are selfconsistent- what enters SDEs is output as well- and the selfconsistency of SDEs 
is more important in QCD then perturbative analyticity, which is not guaranteed in a confining theory at all. In other 
words, what can be obtained by analytical continuation of Euclidean GFs may not agree with the direct Minkowski 
solution. That the other then usually expected singularities appear when making continuation from Euclidean towards 
Minkowski space is a well known fact [1, @. However, in order to be able to correctly check the assumptions, one 
should be able to compare the results based on the continuation of the Euclidean models with the results of the same 
(as possible) model directly solved in the Minkowski space. In this talk I present a direct solution of the SDEs in 
Minkowski space, the first one for QED2+1 , where no analytical assumptions are made. The second model I present 
here is large Nf QCD where some assumptions are necessary in order to make a Minkowski solution possible. 

Before doing so I mention some Minkowski SDEs solution which are based on the assumption of validity of Wick 
rotation. 
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II. SEMI-PERTURB ATIVE METHODS BASED ON INTEGRAL REPRESENTATION 

Assuming perturbative analyticity of Green's function one can impose Khallen-Lehman representation for propa- 
gators 

oo 


with the important property 

oo 

J dk 2 G free {k) = — z7r. (2.2) 

— oo 

Quite independently on the details of the models the momentum SDEs can be turned into a " regular" equations 
for the continuous part of the Lehmann weight p c (M ) (assuming that the full p includes also a single delta function 
corresponding with mass pole of the propagator). Using the Feynman tricks and some standard algebra one can 
always arrive to the dispersion relation for the inverse propagator: 

oo 

G-^Polyn. + J d M 2 k2 ^l.. £ , (2.3) 

T 

where the selfenergy (mass, polarization, etc.) weight function nontrivially depends on p. 

Propagators are complex above the threshold T and the both p and a can be extracted by comparing Re and Im 
parts of G (for a review see Q), the other application of the spectral technique has found its place in the pinch 
technique study of Yang-Mills theories d, [l(| E3 • 

For such Minkowski space calculation the results were always available for small coupling wherein they are quan- 
titatively comparable with the perturbation theory (note, by construction the perturbative contribution is always 
involved in the game) , however for a large enough coupling there is observed some disagreement with the assumption 
or even the solution is not feasible at all Reconstruction of GFs in the spacelike domain, where p 2 < in our 
metric (recall G, p, a are primarily obtained at timelike region (p 2 > T), see @), then using independently p in (|2.ip 
and a in (|2.3[) serves as a simple test of validity of the assumptions. Regarding modeling of QCD GFs the observed 
discrepancy is mild for pinch technique gluon propagator presented in the paper (l3j but becomes a disaster when 
attempting to find a spectral solution for the quark propagator [l2j]. These observations suggest that perturbative 
analyticity is too strong assumption for the QCD quark propagator and very likely for gluon propagator as well (to 
be strict, the position of unusual singularities necessarily affects the analytical structure of the other GFs, since they 
are necessarily related through the SDEs). 



III. QED2+1 

QED2+1 is known [l4j].[l5| to be a simple confining theory which for a small number of flavors posses DCSB. For 
static (heavy) electrons it pos ses a logarithmic confinement if the unquenching effect is not going to spoil the single 
pole of photon propagator [l6tll7l|. QED2+1 is also pedagogical example of numerically soluble theory in Minkowski 
space in ladder approximation of electron SDE [18j | . Here we employ Landau gauge and ^4 = 1 approximation (which 
is exact in Euclidean space, but not obvious or clear in Minkowski space). In addition what I presented in the talk I 
provide a simple derivation for completeness here. 

The full electron propagator reads 

S = S s (p 2 ) i> + S V ( P 2 ) (3.1) 
_ B(p) _ A{p) 

s — 49/ \ 9 7">9 ( \ 1 v 



A 2 (p)p 2 - B 2 (p) ' * A 2 (p)p 2 - B 2 (p) ' 
where in our approximation 



B = m + ^ B , A = 1 ; (3.2) 
E B = -2te 2 J ^L Ss (k 2 )G(q 2 ) , (3.3) 
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where q = k—p ; G = 1/q 2 - We intend to rewrite Eq. (|3 . 3[) into a numerically easily soluble integral equation. In order 
to avoid numerical interpolation we shall use the arguments of propagators as integral variables. Here we perform the 
derivation for the timelike region of p 2 , and without loss of generality we can use the frame where p^ — (p, 0, 0). For 
this purpose we consider a bit more general loop integral. 

/d 3 k 
±^S{k 2 )G{q 2 ){k.q) N , (3.4) 

where q = p — k and where the functions S, G are functions only of variable k 2 and q 2 respectively and clearly N = 
is sufficient in our case. It is obvious that all one loop Lorentz invariant selfenergy contributions can be written in 
terms of considered integral (in any theory). For the space part of the Lorentz three- vector we can use the usual 
spherical coordinates, so the scalar product k.p reads 

k.p = k p - \k\\p\coscj), (3.5) 

and for the measure we have 

/OO /* /•27T 

dk / d\k\\k\ / d<t>. (3.6) 
-oo J JO 

As we choose the argument of S as an integration variable, the integral (|3.4p can be written in the following way 

/ d\k\\k\ * r '£S(k 2 )G(q?)(k.q) N ; (3.7) 

(2tt)^ J 2 \/k 2 + k 2 hi 



q 2 ± =k 2 +p 2 T2p\fk 2 + \k\ 2 ; (3.8) 



(k.q)± =k 2 T P\/k 2 + \k\ 2 ; 

where we have used the explicit form of p^ L . In the next step we split the integral over k 2 to its spacelike and timelike 
subregions and change the order of integration 

/OO />0 />oo />oo />OG 

dk 2 / dk 2 d\k\ + dk 2 d\k\ . (3.9) 

-\k\ 2 J-oo J^W JO JO 

Consequently we make substitution \k\ — » q± in the first (for the "plus" index) and in the second (for the "minus" 
index) term in (|3 . T|) separately, thus getting for / the following expression: 

Jk 2 + P 2 ~4p 



dk 2 / dq%S(k 2 )G{q%y ' W+ 



'• v (V2(q+ + k 2 -p 2 )) N 

k 2 +p 2 -2y/p 2 k 2 ^ + y ' 

dk 2 / dq 2 _S{k 2 )G{q 2 _ y 1 w ~ tLLL- 

-oo Jk 2 +p 2 4p 

"«>r _^_ s y) G y-) (1/2(g - + /- y2))K . (3.io) 

Jk 2 +p 2 +2y?p 2 k 2 4p 

In what follows we use a more compact notation by relabeling </ 2 — > q 2 and q 2 . — > q 2 , further we can rewrite the 
boundaries in a fully equivalent manner 

m = r * 2 r ^^( fc 2 )G(g 2 ) M^i^i^ 



— oo J — oo 

oo 2 ,fcW+2V?F (1/%2 + fc2 _ p2))J V 



dfc 2 / rfg 2 ff(fc 2 )G( g 2 ) v ' v * ' Z-U— (3.11) 

k 2 +p 2 -2jp 2 k 2 —4p 



which is the most general expression whether G is known or not. 
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G is the free boson propagator in our approximation for which < G >= 0, so the first term vanishes (we did not 
consider ie prescription here, as it has no important effect in this case) and the q 2 integration can be performed 
analytically leading to the final expression for (|3.3p 

ie2 f°° jJo ^v. I ^ — ^1 



B(p') = m +— dk z S s {k z )ln^ ^. (3.12) 

4?H J (k + p) 

Recall that the same has been derived in [l8[ by using a hyperbolic coordinates. 

We assume the propagator functions are complex almost everywhere (i.e. there are no perturbativc thresholds) 

B(p) _ RbCi+T b c 2 T B ci~R B C2 f , qiq\ 

s ~ A 2 (p)p 2 - B 2 (p) ~ cj + 4 +l cl + cl ; ( ' 

Rac\ + T A c 2 T A ci - RaC2 

v 2 I 2 ' 2 I 2 

c l i c 2 C l "T" c 2 

d = (^-r^) P 2 -(i?|-r|); 
c 2 = (i? A rV - J2 B r|) ; 

where R, T are Re an Im parts of proper GFs, A = Ra + i^A ', B = R B + iF B . 

In this way we get two coupled integral equation for R Bl T B which can be solved by the method of iterations. These 
integral equations are regular in S unless T B does not vanish. 



IV. QED2+1 VERSUS QED3 

At this point it is quite interesting to compare formally obtained expression for Minkowski B with its Euclidean 
counterpartner. Stress the both variables k,p in (|3.12p are timelike fourmomenta since the timelike solution decouple 
from the spacelike one. Furthermore as we will see S s given by Eq. (|3 . 1 3|) is entirely complex function in the timelike 
Minkowski subspace. 

The conventionally written Euclidean partner (note p 2 E = —p 2 > ) is a smooth regular function 

SsiP%) = ^b|K + Wl) ■ 
When using the same approximation it satisfies Wick rotated Euclidean ladder SDE: 

oo 

B E {p%) =m-^J dk 2 E S s {p 2 E )ln^-^ (4.2) 
o 

where B B is numerically known. It is real, and it stays real even if one allows for a complex valued propagator 
(imaginary part is not generated), B is known to be nontrivial for any coupling e which includes classical chiral 
symmetry (to = 0) case. 

Due to the very simple structure of SDEs we can find recipe how to get the first from the second one, however one 
cannot say that the first is the analytical continuation of the second. Very formally, the Euclidean ladder SDE can be 
get by using the contour shown in Fig.l. with simultaneously continuing external variable p to the spacelike region 
as well (explicitly pe —> —ipM- 

Stress here, this continuation is performed without further deforming the contour in Fig. [T]and considering contri- 
butions from complex branch points observed in In other words, we know the prescription for the transformation 
of equation between two different spaces- Minkowski and Euclidean-, however the timelike Minkowski solution is not 
an analytical continuation of the Euclidean one. The numerical results on Minkowski QED2+1 have been obtained 
first time in the paper [l8j and here it is presented in figure [2] and [3] Apparently, the mass function B would be 
non-holomorphic in the beginning p 2 — if one tries to interpret the Euclidean solution as an Minkowski spacelike 
solution. Therefore the Wick rotation is not valid here. In the limit e 2 >> to we get chiral symmetry breaking phase 
<t>M D csB = 88°. 

Until now a direct Minkowski solution of B for spacelike argument is unknown to us (attempting to put similar 
ideas in a game, e.g. using the arguments of GFs as the integration variables, we would certainly factorize some 
singularities that remain isolated in different subregions of Minkowski space). However to deal with this problem, one 
can employ correct analytical continuation of the obtained timelike solution (so the strategy is completely opposite 
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FIG. 1: Graphical explanation of Minkowski -£ Euclidean continuation for ladder QED2+1. The figure shows the complex 
plane of the integral variable \fW used in the 1-dim integral SDEs in Minkowski and Euclidean spaces. The recipe that gives 
the correct continuation is based on the integration along the drawn contour such that we take J Imk — * J Rt , k we going from the 
first to the second space, albeit Cauchy lemma for such contour cannot be used, since the function is not holomorphic in the 
interior of the shown curve. The arc does not contribute as radius goes to infinity because of Jordan lemma. 



Magnitude of dynamical mass in QED2+1, m=l 
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FIG. 2: Magnitude of electron dynamical mass function 

to the standard "text book" procedure, here we get the timelike solution at first and then we can construct the 
solution at spacelike region). In given approximation the analytical continuation is a very straightforward if one takes 
In —>■ Ln. The resulting equations reads 

B(-p 2 ) =m+ [ dk 2 S s {k)arcsin J*? (4.3) 

v 

where S s in the rhs. of Eq. (|4.3[) is defined in timelike region. The numerical integration should also be straightforward 
and will be done elsewhere. 

The second very important observation is that we get no real pole in the electron propagator. Actually B is 
complex for all considered (real) p 2 and if bare mass m is not very large it has relatively large non-zero imaginary 
part preventing thus the presence of a real pole. No real pole implies no free particle solution, which is the simple 
way how confinement of QED2+1 electron is reflected by behaviour of the GFs. 



Dynamical mass phase of QED2+1 electron 
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FIG. 3: Phase of dynamical electron mass M — \M\e'' 



V. DCSB IN LARGE N f QCD 

Up to date there is no known direct nonperturbative Minkowski space study of a pure Yang-Mills theory, e.g. N c = 3 
pure gluodynamics. There are attempts in the literature that at least questioned the analytical structure of GFs in 
QCD. Serious study of analytical continuation of the Euclidean quark propagator performed in Q shows branch 
points in the complex plane of fc 2 . The branch points are located in the interior and as well as in the exterior of Wick 
rotation contour, however the positions and number of singularities found depends on the details of the interaction. 
Recall also that the character of singularities is unknown in general, very likely they are not simple isolated poles. The 
truncation of SDEs is a necessary approximation in any case and discussion of the effect of a truncation is plausible, 
however it seems quite natural that a strong enough interaction likely leads to the absence of a real branch point 
in the quark propagator. Also the naive numerical fits and indirect estimates based on the behaviour of Schwinger 
functions made for instance in 0] signaled that the quark propagator has a more complicated analytical structure. 
Note also that that the fit used in Q is in contradiction with the Wick rotation. As in the previous case I am going to 
discuss the fermion SDE, but now in 3+1 Minkowski space. Using a certain simplification the model has a confining 
phase characterized by absence of a real fermion propagator pole. In this respect we are confirming proposals by 
Fukuda, Kugo but reaching quantitatively very different result since avoiding using Euclidean space calculation 
at intermediate step. 

In order to get some first solution we consider large Nf QCD. Here the number of massless quarks is larger then in 
the nature, decreasing considerably low q 2 running, but limited from the up in order to preserve asymptotic freedom. 
Here we assume that Euclidean theory is at least a good guide for the estimate of the magnitude of an effective 
running coupling. So for instance Nf — 7, 8 is a reasonable estimate for N c = 3. For this number we still should get 
chiral symmetry breaking. Typically in this and similar models the form factor and related effective coupling do not 
run but mildly change in the low scale (defined by the low energy mass of the quarks m « A) and vanishes at high 
scale in accordance with asymptotic freedom. We basically use this scenario here we consider the same model as in 
the paper [2l|. It is also notable, these Technicolor (TC) like model remains viable scenario for electroweak symmetry 
breaking without inclusion of the scalar Higgs field [22j • 

In our Minkowski study we model the effective running charge as follows 

for \q 2 \ < A 2 : a TC (q 2 ) = a* ; 
for | g 2 | > A 2 a T c(<7 2 H0 , 

where constant a* is large enough to get DCSB. This coupling enters the ladder approximated SDE which in Landau 
gauge and A = 1 approximation reads 

M( P 2 ) = m + icJ^S s (k)^ 1 - (5.D 



The general structure of the quark propagator is given again by (|3 . L3|) and the dynamical mass function is simply 
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M = B. C is the constant stemming from Casimir of given representation of Non Abelian gauge group and the 
prefactor from the Lorentz and Dirac algebras. 

Minkowski momentum space integral which appears in our SDE is just an one more spacelike dimension extended 
analogue of / considered for QED2+1, here it explicitly reads 



G(k 2 )S(q 2 ) 



(5.2) 



where G = arc i 9 m our case. 

Lorentz invariance dictates that B(p 2 ) is the function of onlyp 2 for any choice of p M , so to get a solution it seems 
to be advantageous to use some simple choice of configuration, for spacelike p M = (0, 0,0, p) while for timclikc p 2 
p» = (p, 0,0,0). 

In the usual Euclidean studies it is desired to use arguments of GFs as an integration variables. Unfortunately for 
a spacelike p it leads to factorization of an awfully divergent or singular integrals which are completely independent 
on the behaviour of GFs inside the integrals . Actually one can arrive to the following singular integral 



for p 2 < 0, 
I[G,S;p)~ 



1 



1 



dz 



( 



\ 



oosgn(z) 



dk 2 



z 3 k 2 + a/ z 2 
(q 2 -k 2 



K = 



k 2 +p 2 

KG(k 2 )S(q 2 ) 

p 2 ) 2 . _ (q 2 k 2 



oosgn(z) 



dk 2 



2\2 



V) 



so the problem is very badly defined. It is almost redundant to say that in perturbation theory this is just the Wick 
rotation with its ie prescription which regularize unwanted singularities. We do not know a different Minkowski space 
rcgularization which preserves all required symmetry (Lorentz and gauge if required) for spacelike p 2 , however for 
positive p 2 one can arrive to a more optimistic formula, which when used in our SDE (|5.1[) gives the following integral 
equation: 



M(p 2 ) = 



iC 



16p 2 7T 3 



dk 1 



(fe-p) 2 



dq 2 V / A(q 2 ,k 2 ,p 2 ) + / dq 2 y / A{q 2 ,k 2 ,p 2 ) 



(fc+p) 2 



/U POO 
dk 2 / dq 2 v / A(q 2 ,-k 2 ,p 2 ) 
-OO J —OO 



a T c{q 2 ) 



(5.3) 



A(a, b, c) = a 2 + b 2 + cf - 2(ab + ac + be) 



Not similarly to QED2+1 the equation (|5.3p remains coupled with unknown S s for the spacelike arguments. It 
requires some further approximations which together with the detail of the derivation of Eq. 15.31 will be published in 
the shortly coming paper [2fJ. 

The present numerical solution has been obtained for the effective coupling a* = 1.5 and the constant prefactor 
C has been absorbed into it. Until now we did not perform an ultimate study of large parameter space, but we can 
conclude that we have never observed nontrivial solution for exactly zero bare mass m = 0. Since the value can be 
arbitrarily small we were slowly decreasing m during the iteration process to a desired value. The bare mass we used in 
our calculation is presented in the figure [5] All results are scaled with respect to A which is quite natural scale (recall 
experimental Aq C£) = 450MeU while expected ^-Technicolor = 1 — lOTeU. The solution we get is completely stable 
and in principle we can reach it with an arbitrary high accuracy. The most apparent fact is the infrared behaviour, 
the mass simply goes to infinity as l/p there, in other words constituent (techni)quark mass (even as absolute value) 
is not well defined quantity since lim p 2^ Q ReM{p 2 ), ImM(p ) — * oo. At the time being we do not know whether this 
is a pathological feature of the Minkowski metric in used, or it it is a physically acceptable scenario for strong coupling 
3+1 dimensional theories. This remains to be a subject of future studies and confirmations (based for instance on the 
bound state solutions). 
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FIG. 4: Left panel shows real and imaginary parts of the dynamical mass separately, the right panel display the same as 
previous figure but in log-log scale for the both quantity, mass is rescaled by A, 




FIG. 5: On the left panel we plotted the magnitude <j> (in rad) of dynamical quark mass in 3+1 Minkowski dimension, the 
phase 4> is on the right, note M = ||Af||e i0 . 



Further, like in QED2+1, there is no physical pole mass, the mass function is complex preventing the existence of 
a real pole of the propagator. Following the LSZ-reduction formula for S-matrix, the matrix elements between quark 
states is zero, since the limit lim p 2^ m 2(p 2 — m 2 )S(p 2 ) = for any considered real m 2 , or in better words: free quarks 
do not exist. 

To conclude in one sentence, there are many issues to be clarified in future Minkowski space nonperturbative studies, 



but the clue to be followed already exists and the perspectives are slowly opening. 
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